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Abstract 


The application of complex networks of elementary chemical reactions to the gas phase of 
burning energetic materials has increased markedly over the last decade. The exquisite 
complexity of these gas-phase reaction networks, coupled with available high-rigor treatments 
of transport, is not matched by an equivalent level of sophistication in descriptions of the 
condensed-phase and interfacial phenomena. Owing to the vastly more complicated, many-body 
nature of the condensed phase, this condition is not likely to be relieved soon. In response to 
these difficulties, a new semi-empirical approach to burning-rate calculation has been developed 
and applied to frozen ozone, cyclotrimethylenetrinitramine (RDX), and nitroglycerine. The new 
approach hypothesizes a single overall reaction linking the unreacted material to the net products 
of condensed-phase decomposition and characterizes their rate of formation according to an 
empirically derived pyrolysis law. These condensed-phase products become the reactants for the 
gas phase, which are treated in full elementary-reaction detail. Using this new semi-empirical 
model, a methodology for computing the relative effects of several additives on the burning rate 
of nitroglycerine is developed and demonstrated. Hopefully this approach will enable more rapid 
progress in modeling multi-ingredient propellants than did previous approaches attempting to 
model the condensed-phase processes in detail. 
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1. Introduction 


The last 10 years have seen an important paradigm shift in the modeling of energetic material 
combustion. Models using chemically nonspecific reaction formalisms have given way to models 
that promise to explain the detailed chemical nature of combustion through their explicit elementary- 
reaction mechanisms. The seminal prototypes for the chemically nonspecific models were 
developed during WWII by Parr and Crawford [1] and Rice and Ginell [2]. Typically these models 
considered a gasification reaction at the surface of the condensed phase and a single overall reaction 
in the gas phase; burning rates were computed by imposing the energy and species conservation 
equations in each phase and at the burning surface. A conceptual review of many of these models 
was given by Miller [3], and a review of their application to composite propellants was provided by 
Ramohalli [4]. These early models elucidate the interplay between the chemical and physical 
phenomena that produce the observed patterns of burning rate dependence on pressure and initial 
temperature and provided guidance in the choice of some physical features of composite-propellant 
formulation, such as oxidizer particle size. By their nature, however, these models cannot give 
guidance for chemical formulation and the effects of chemical additives on performance. This role 
had to await the development of chemically specific descriptions of the combustion. Models with 
explicit chemical mechanisms now dominate new development activities and are the subject of this 
article. We will confine ourselves in this report to steady-state combustion, although some recent 
work [5, 6] has addressed the role of elementary reactions in ignition. 

The transition from overall reactions to detailed elementary reactions in the modeling of 
energetic material combustion did not occur suddenly and, in fact, is not yet complete. Thus far, the 
use of elementary reactions has been limited to the gas phase; our fundamental understanding of the 
condensed phase processes is not allowing more sophistication at present. In the interim, 
condensed phase reactions have been treated as one or more overall reactions. Reaction mechanisms 
for the gas phase were first developed apart from any attempt to compute the burning rate. As early 
as 1965, Softer [7] assembled a mechanism consisting of 12 species reacting by 17 reversible 
reactions to describe the secondary reaction zone of double-base propellants. He numerically 
integrated this mechanism, neglecting transport processes, and found order-of-magnitude agreement 


1 



with experimentally determined induction times for the visible flame. In the 
early-to-mid 1980s, Ermolin et al. [8] worked out elementary-reaction mechanisms for ammonium 
perchlorate and cyclotrimethylenetrinitramine (RDX) [9], guided by mass spectrometric 
measurements in the flames of these condensed systems. 

Attempts to compute the burning rate of an energetic material using realistic chemistry probably 
began with the Guirao and Williams’ [10] treatment of ammonium perchlorate in 1971. Though the 
chemical mechanism consisted of 14 species and 10 irreversible reactions, a small problem by 
today’s standards, computational tools at that time did not permit the simultaneous solution of the 
kinetics and transport. A simplified simultaneous solution of kinetics and transport for a larger 
mechanism of 19 species and 60 reversible reactions was reported by Hatch [11] in 1986 for 
nitroglycerine. A significant advance in the computational tools available for work of this type 
occurred in 1985 with the publication and availability of the PREMIX code [12], a user-friendly, 
well-documented, one-dimensional (1-D), premixed, laminar flame code. This code and its later 
improvements provide a high level of rigor to the description of gas-phase transport, allowing for 
such subtle effects as thermal diffusion and multicomponent transport to be conveniently treated. 
Moreover, the description of thermodynamic functions, transport parameters, reactions, and reaction 
rates in this code is also very general, allowing a wide scope of kinetic representations and 
convenient updating. Melius [13] was the first to utilize the PREMIX code in a calculation of 
energetic-material burning rate, addressing the case of RDX. In the United States, as a result of 
several Office of Naval Research and U.S. Army Research Office workshops starting in about 1987, 
a concerted effort was directed at understanding the physics and chemistry of RDX as a prototypical 
system. This focus resulted in a number of further refinements of RDX modeling [14-17]. 
Extensions of these general methods have been made to include cyclotetramethylenetetranitramine 
(HMX) [18, 19] and glycidyl azide polymer (GAP) [20]. In order to study the physical aspects of 
three-phase combustion. Miller [21] described a detailed chemical model applied to the combustion 
of frozen ozone as a prototypical energetic material for which very reliable kinetic and 
thermophysical data exist. 
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This report focuses on models that compute the burning rate based on detailed elementary 
reaction mechanisms in the gas phase. However, it is worth mentioning a few models that are 
chemically specific to some degree and utilize overall reactions to describe the gas phase. An early 
example of this type of model comes from BenReuven et al. [22] and BenReuven [23]. Their model 
treated the combustion of RDX using one overall reaction in the liquid phase, evaporation as the 
surface-gasification mechanism, and two overall reactions in the gas phase. Another example was 
advanced by Bizot and Beckstead [24], who considered double-base propellant combustion using 
three overall reactions in the condensed phase and two overall reactions in the gas phase. A recent 
model attempting to simplify the very complex chemistry involved in double-base propellants is that 
by Song and Yang [25]. Like the Bizot and Beckstead model. Song and Yang lump species into 
categories such as aldehydes and oxidizers. Yet another approach in this vein was taken by Li and 
Williams [26], who argued for the dissociation of HONO as the rate limiting step in the primary 
gas-phase mechanism of RDX and then applied an asymptotic analysis to compute the burning rate 
based on this one reaction. The relative simplicity of these models is appealing, but it will take time 
to determine the extent to which generality and even understanding is sacrificed. 

2. Conceptual Framework 


As previously stated, our primary interest here is in that class of models which strive to describe 
the gas-phase reactions as true elementary reactions, the rate coefficients for which are determined 
from independent kinetics experiments, reaction-rate theory, or thermochemical estimation methods. 
This restriction narrows the number of models to only four [14,16,17,21]. In this report, we will 
describe a fifth approach. In developing this new approach, we will apply it to RDX. It will be 
helpful, therefore, to review some of the details of the other RDX models [14,16,17] (i.e., those in 
our restricted subset of model types). 

All three models [ 14,16,17] defer to the experimental work of Brill et al. [27] for a description 
of the condensed-phase reactions in RDX. Brill et al. [27] found that RDX decomposes by two 
competing channels, R1 and R2 (see Table 1). Unfortunately, his experiment was unable to 
distinguish between several possible R2 paths, R2a, R2b, and R2c. Both Prasad et al. [16] 
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Table 1. Liquid-Phase Reactions in Deflagrating RDX Adopted by Various Models 


Reaction 

Number 

Reaction 

Heat of 
Reaction 
at 570 K a 
(kcal/mole) 

Model 

R1 

RDX -> 3 CH 2 0 + 3 N a O 

-47 

All [14,16,17] 

R2a 

RDX -* 3 H 2 CN + 3N0 2 

179 

P&S [16], 

D&B [17] 

R2b 

RDX -» 3 HCN + 3 HONO 

19 

— 

R2c 

RDX -» 3 HCN + 3 N0 2 + 3 H 

257 

— 

R2d 

RDX -» 3 HCN + 3/2 NO + 3/2 NCh + 3/2 H 2 0 

34 

L&Y [14] 


a This is the completed surface temperature at 1 atm. 


and Davidson and Beckstead [17] elected to use R2a for the second path. Liau and Yang [14] 
suggested and used in their model another interpretation of Brill’s measurements, R2d. The table 
shows that the degree of endothermicity among the different R2 paths varies dramatically. Unless 
condensed-phase reactions are unimportant to the deflagration of RDX, one would expect such 
differences to have a significant effect on the computed burning rate. Because the three models 
being considered differ in many details other than the condensed-phase reactions just discussed and 
the burning rates computed by each of the models agree quite well with experimental burning rates, 
it is impossible to see the effects of the different liquid-phase reactions in isolation. Nonetheless, it 
may be significant that the amounts of RDX decomposition reported to occur in the condensed phase 
vary among these models from 40% [16] to 25% [17] to a “limited” amount [14]. 

Even though the computational power now exists to treat condensed-phase reactions, our current 
knowledge of these reactions (even at the level of overall descriptions, let alone an elementary- 
reaction description) is very limited. This uncertainty is not confined simply to the rate coefficients, 
but extends to the very identities of both the reactants and products involved in these reactions. We 
believe that these difficulties will not soon be resolved, and it is with this prospect in mind that we 
propose the following new and more practical approach. 
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A key difficulty with chemically specific models is determining the identity and information 
rates of mole fractions of chemical species which first emerge from the condensed phase into the gas 
phase. We will show that given this knowledge and a reaction mechanism for the gas phase, the 
burning rate can be computed. One additional requirement is an expression relating the burning rate 
to the surface temperature; this will be discussed in this report. Though there is speculation required 
in identifying the chemical species emerging from the condensed phase, it is much reduced from 
having to speculate on the reactions occurring in the condensed phase and the rate coefficients of 
those reactions. Techniques to measure or calculate those condensed-phase processes are not likely 
to be realized soon; thus, the semi-empirical approach described here can be thought of as a 
practical, interim strategy pending the development of definitive tools for investigating the 
condensed phase. With these enabling simplifications, the model can be readily adapted to many 
worthwhile purposes, including experimental testing of gas-phase reaction mechanisms, computing 
the burning rate of multi-ingredient propellant formulations, and predicting the effect of chemical 
additives on the burning rate. 

3. Mathematical Framework 

We consider the 1-D, steady-state deflagration of a condensed substance, which is oriented in a 
coordinate system moving with the linear regression velocity such that the unreacted material at an 
initial temperature of To extends to x = - », the regressing surface is always at x = 0 with 
temperature T s , and the final gaseous combustion products at temperature Tf are found at x = + «>. 
The mass conservation equation can be integrated simply to give 


m = pu = p s r , (1) 

where m is the mass burning rate or mass flux, p and u are the mass density and mixture mass 
velocity at any point in the interval, respectively, p s is the mass density of the unreacted material at 
its initial temperature, and r is the linear regression rate of the surface. The continuity equation for 
each chemical species i is 
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(2) 


-f (Wi) + ~ m 

dx ax 


0 . 


Here F; is the mass fraction of the zth species, V, is the diffusion velocity of the ith species, is the 
molar production rate by reactions of the zth species per unit volume, and Wi is the molecular weight 
of the zth species. Finally, the equation of energy conservation for this system is 


Mi*' 

dx dx 




N 


N 


me, - Ipyy.c,, - 

uX ; i 


0 , 


( 3 ) 


where X is the mixture thermal conductivity,^ is the mixture specific heat, c/>, is the specific heat 
of the zth specis, hi is the enthalpy of the zth species, and N is the total number of species being 
considered. 


These equations are to be solved in the condensed and gas phases and their solutions matched 
through appropriate boundary conditions at the surface. In the condensed phase, the energy 
boundary conditions are 


T = T 0 at x = — oo 

T = T x at jc = 0, (4) 

and the species boundary conditions are 

Y i = Y~ at x = -oo 

Y i = Y~* at x = —0. (5) 

In the gas phase, the energy-equation boundary conditions are 
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T =T S at * = +0 


dT n 

— = 0 at x = +oo, 
dx 


( 6 ) 


and the gas-phase boundary condition on each species is 


Y ; = 1?° at x = +0 


dY ; 


l _ 


dx 


= 0 at x = +oo. 


(7) 


At the surface between these two regions, the continuity of species and energy flux is guaranteed by 
the following two boundary conditions. Continuity of species flux at the surface is expressed by 


nX? + = ml ? 0 + pT.A? 0 , 


( 8 ) 


where the superscripts -0 and +0 and subscripts c and g refer to the condensed and gas sides of the 
surface, respectively. Continuity of the energy flux across the surface boundary is expressed by 


( / fr \"° n n 

—J + mXirv + p, trrv.v = 


-A. 


dx 

dT 


J 

fsr V° 


K dx j 


+ m^r, a h“ + P t *“*,*“■ 


(9) 


These general equations are the starting point for our simplifications. In the gas phase, these 
conservation equations are solved numerically using the PREMIX [12] code as a subroutine; thus, 
the gas phase is treated with full rigor in our model. In the condensed phase, we assume that 
(1) there are no in-depth chemical reactions, and (2) molecular diffusion is negligible. With these 
assumptions, the energy equation in the condensed phase reduces to 
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( 10 ) 


dx 


dT') 
X — 
dx) 


- m- 


dx 


f N \ 

E Y A 

i J 


= 0. 


Integrating this equation from x = - oo tox= : -0 results in 


K 


d T 
dx 


\-° N 


= m^Y: 0 h; 0 - mYjrK 


Substituting equation 11 and equation 8 into equation 9 yields 


X. 


d T 

\dx j 


= mi(w - 


(id 


( 12 ) 


Equation 12 is the form of the energy-flux boundary condition used in the model developed in 
this report. Notice that one does not need to know the enthalpies of the species on the condensed 
phase side of the surface at T s , only their mass fractions there. This is important since heats of 
desorption of the nascent gas-phase species may be difficult to estimate. We are unaware of 
previously published material on this form of the energy-flux boundary condition, perhaps because 
of the novelty of the present context of a single overall condensed-phase reaction being coupled to a 
full elementary-reaction description of the gas phase. This equation enables our treatment. Notice 
also that neither the thermal conductivity, the mass density, nor the specific heat of the condensed 
phase at temperatures near the surface are required; the values of all these quantities are very 
uncertain. Only the specific heat and mass density of the unreacted material over the range of initial 
temperatures are required, at least for purposes of computing the burning rate. Of course, if the 
temperature profile through the condensed phase is needed, then a further integration using values 
for these parameters over the full condensed-phase temperature range will be required. 
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In order to compute a burning rate, the following information must be supplied; the enthalpy and 
mass density of the unreacted energetic material over the range of initial temperatures, a set of 
products (species and mole fractions) of the condensed-phase decomposition, a rate of appearance of 
these decomposition products as a function of surface temperature, and an elementary reaction 
mechanism through which these decomposition products react in the gas phase. The condensed- 
phase decomposition products become the initial or nascent gas-phase reactants. 

As pointed out in our conceptual framework discussion, the identities of the condensed-phase 
reaction products are not known completely, even for the relatively well-studied case of RDX. For 
most other energetic materials, the situation is worse still. The only clear constraint we have in 
constructing this species set is elemental balance. Therefore, from the universe of possible product 
sets (i.e., those leading to a balanced overall reaction), plausible products must be selected based on 
either experimental knowledge of these product identities and concentrations (if available) or 
theoretical consideration of likely reaction paths. In general, multiplicity of such plausible product 
sets may occur. In this case, that which best reproduces the experimental burning rate should be 
chosen. Or, an average of the product sets can be taken based on the assumption of equal a priori 
probabilities of occurrence. In any case, the decomposition set chosen can be improved as new 
experimental data or theoretical insights become available. 

The final type of information required to compute a burning rate is the rate of appearance of the 
condensed-phase decomposition products. For this, we appeal to the very old notion of an 
Arrhenius-like expression [2] relating the burning rate to the surface temperature, such as 


m = A s e 


-EJRT, 


(13) 


Such an expression has often been termed a “pyrolysis law” in the propellant-combustion literature. 
There are other forms of this expression in use [28,29], but the above form preserves the dominant 
functionality without additional (probably uncertain) data. In fact, Zenin [28] prefers equation 13 to 
represent his experimental determination of surface temperatures vs. burning rates for a range of 
double-base propellants. Zenin found that this form provides a universal representation for 
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double-base propellants of different ingredient proportions (see Figure 1). It is this universality for a 
given class of propellants that impels us to adopt it here as a formalism for the rate of appearance of 
the condensed-phase decomposition products. In the case of double-base propellants, the 
mechanism for surface regression is likely to be reactive in nature because of the polymeric character 
of the nitrocellulose molecules. Figure 2 shows that the pyrolysis law in equation 13 also proves 
appropriate to a purely evaporative surface regression mechanism. Here, the surface temperature for 
the three-phase deflagration of frozen ozone was computed based on nonequilibrium evaporation 
driven by heat feedback from gas-phase reactions [21]. The steady deflagration of RDX is also 
believed to be predominately the result of evaporation [14]. Figure 3 shows a best-fit pyrolysis law 
(in the least-squares sense) for RDX using the surface-temperature and buming-rate data of Zenin 
[28]. Again, the fit is quite reasonable. 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 

10 3 x1/T s (K) 


Figure 1. Demonstration of the Universality of the Pyrolysis Law for Double-Base Propellants 
Based on the Work of Zenin [28]. 
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4 5 6 7 


10 3 x1/T s (K) 

Figure 2. Quality of Fit of a Pyrolysis Law to an Evaporative Regression Mechanism 
Calculated From a First Principles Treatment of Frozen Ozone Deflagration [21]. 



1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 

10 3 x1/T s (K) 

Figure 3. Pyrolysis Law for RDX Derived From Zenin’s [28] Experimental Data. 
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4. Example: Frozen Ozone 


Our first example of the application of this new model is that of the self-sustained deflagration of 
frozen ozone. Using the pyrolysis law previously determined from a first-principles model [21] to 
test for consistency, the burning rate is computed by the new model and compared with the burning 
rate computed by the first-principles model. Such a test is not as trivial as it might seem. The 
calculational details involving the evaporative surface-regression mechanism of the first-principles 
model is considerably more complicated than that involved in the new model. However, the 
evaporation process is characterized by an energy barrier (heat of vaporization), and the assumption 
of local thermodynamic equilibrium introduces a Boltzmann factor, this suggests that the process 
may be controlled by an Arrhenius-like term. On the other hand, the evaporation-mechanism 
equations [21] do not lead to such a simple expression, though an Arrhenius-like dependence on 
temperature evidently does have an implicit relevance given the excellent fit of the pyrolysis law 
observed in Figure 2. 

Since any condensed-phase reactions in liquid ozone are too slow to influence the deflagration 
rate, the mechanism of surface regression is the following: evaporation driven by heat, released in 
gas-phase reactions, and conducted back to the surface. In this situation, the nascent gas-phase 
species is evaporated ozone (i.e., = 1). Coupling this with the fitted pyrolysis law just described 

and the reaction mechanism consisting of three reversible reactions [21], we have computed the 
frozen ozone burning rate by our new model. It will be compared to the rate computed by the first- 
principles model [21]. Figure 4 shows that the agreement is excellent. The first-principles model in 
this case includes thermal diffusion and multicomponent transport, and the numerical grid spacing is 
fine enough that the burning rates are computed to an accuracy within a few tenths of a percent. 
Also, the central differencing option was selected. Thus, the success both of the pyrolysis law and 
the new model is based on a high level of computational rigor. The dominant chemical steps 
occurring in the gas phase for the ozone case are thoroughly discussed in Miller [21]. 
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As previously discussed, the condensed-phase decomposition paths for RDX are uncertain. 
Evaporation was determined to be the overwhelmingly dominant mechanism by Melius [13], Liau 
and Yang [14], and Liau [15], and at least predominant by Davidson and Beckstead [17] and Prasad 
et al. [16]. We therefore examine the consequences of assuming in our model that the surface- 
regression mechanism is purely evaporative. Thus, by analogy to the ozone case above, vapor-phase 
RDX is assumed to be the sole species emerging from the surface; mathematically this is expressed 
as Y^x = 1 • The pyrolysis law is obtained by a linearized least-squares fit to the experimental data 
of Zenin [28] and is illustrated in Figure 3. 

The enthalpy of the solid at the initial temperature is determined as follows: the specific heat of 
solid RDX is obtained by a least-squares fit of a linear temperature function to the data of 
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Shoemaker et aL [30] and Miller [31]. The fit is shown in Figure 5 and its parameters are in Table 2. 
This functional form of the specific heat is integrated over temperature to determine the enthalpy of 
solid RDX using the value of the heat of formation at 298 K to obtain the integration constant. This 
treatment differs from that of Liau and Yang [14] and Liau [15], in that their model assumes a 
constant value of specific heat. 
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Figure 5. Linear Least-Squares Fit of the Specific-Heat Data of Shoemaker et al. [30] and 
Miller [31] for Solid RDX Over a Range of Initial Temperatures. 


Table 2. Values of Condensed-Phase Parameters Used in RDX Calculations 


Value 

Description 

As = 6.134 x 10 5 g/cm 2 -s 

Exponential prefactor in pyrolysis law 
(see text) 

Es- 18,539 cal/mol 

Activation energy in Pyrolysis law 
(see text) 

ps= 1.66 G/CM 3 

Mass density of pressed solid used in 
experiments [28] 

Wrdx = 222.118 g/mol 

RDX molecular weight 

AHf K = 14,590 cal/mol 

RDX heat of formation at 298 K 

c“ Kd =0.03604+ (7.105 x 10‘ 4 )T cal/g-K 

Specific heat of solid RDX over temperature 
range 260-400 K (see text) 
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The reaction mechanism (Yetter’s second model [16]), thermodynamic data, and transport data 
for this case is taken from the work of Liau [15]. The burning rate (Figure 6), temperature 
(Figure 7), and species profiles (Figure 8) computed by the present model at an initial temperature of 
293 K agrees very closely to that computed by the Liau model. (Dr. Y.-C. Liau kindly made his 
code available to us. All calculations in this article attributed to the Liau model [15] were made 
using the code version dated 13 February 1997.) We emphasize that the gas-phase reaction 
mechanisms, species thermodynamics, and transport parameters in both models are identical; only 
the treatment of the condensed phase and surface-gasification mechanism are different. The model 
burning rates shown in Figure 6 also compare favorably with the experimental data of Miller [32], 
Zenin [28], Ulas et al. [33], Homan et al [34], and Atwood et al [35]. 



Pressure (atm) 


Figure 6. Comparison of Burning Rates for RDX Computed by the New Model to Those 
Computed by That of Liau [15] and to Experimental Data. (See Text for 
Experimental Data References.) 
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Distance from Surface (microns) 

Figure 7. Temperature Profiles for RDX at 1 atm and 293 K, as Computed by the New 
Model, Compared With That of Liau [15]. 

In the RDX case, we have sought to minimize nongermane computational differences between 
the new semi-empirical model and that of Liau and Yang [ 14 ]. Their choice of the 
windward differencing option, numerical grid options, and suppression of the thermal diffusion and 
multicomponent transport options in PREMIX were selected for our calculations as well 

The dominant reaction paths for RDX flames have been discussed previously [ 13 ], and while the 
reaction mechanism has been expanded and refined since that study, we believe that the conclusions 
have not changed significantly. Hence, no new dominant-path analysis for RDX is given here. 
However, a brief explanation is now given for the slight rise in temperature above the adiabatic 
equilibrium value late in the flame zone, as shown in Figure 7 . This behavior is surprising compared 
to 02/fuel combustion in which the temperature rise is monotonic; however, it is quite typical of 
flames oxidized by N2O and NO2. In O2 flames, there is typically an overshoot of H, O, and OH 
concentrations late in the flame zone, above their eventual equilibrium values. In hydrocarbon/02 
flames, this overshoot is caused by the following radical chain-branching reactions: 
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Mole Fraction 



Distance from Surface (microns) 

Figure 8. Species Profiles for RDX at 1 atm and 293 K, as Computed by the Present Model 
(Solid Lines), Compared With That of Liau [15] (Dashed Lines). 
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H + Oz -> OH + O 


and 

CH 3 + O 2 CH 3 O + O. 

As this excess of radicals recombine in the bumt-gas region, heat is released and the temperature 
continues to rise to the adiabatic equilibrium value. For N 2 O- and NCVoxidized flames, these 
chain branching reactions are absent. Towards the end of the flame zone, the radical concentrations 
are usually somewhat below their eventual equilibrium values. To reach the final equilibrium 
concentrations of radicals, stable species slowly decompose in the burnt gas region, absorbing heat 
in the process, and lowering the final temperature to the equilibrium value. For the case in 
Figure 7, our calculations show that H, O, and OH concentrations rise by factors of 2 to 5 between 
800 and 2,000 |im. Over this same region, the temperature drops by about 100 K. A similar result 
was observed in computations for a low-pressure H 2 /N 20 /Ar flame [36]. 

6. Example: Nitroglycerine 


Relatively little information is available on the thermal decomposition pathways of 
nitroglycerine (NG). We examine it here in the context of the new semi-empirical burning rate 
model because it represents a case of intermediate chemical complexity, it continues to play an 
important role in gun propellants, and it is a substance for which experimental burning rate data 
exists. Vital also is the fact that we could assemble a gas-phase reaction mechanism with relatively 
few changes from that previously developed for the dark zone of nitrate-ester and nitramine 
propellants [37,38], This was accomplished by adding only one additional species, CH 2 O, and its 
associated reactions. The resulting mechanism (labeled DB11 here) is given in the Appendix and 
consists of 35 species and 178 reversible reactions. Other parameters are given in Table 3. 

For the pyrolysis law, we simply use the universal double-base law determined by Zenin [28] 
and shown in Figure 1. Double-base propellants, of course, consist primarily of different proportions 
of NG and nitrocellulose. Though NG is therefore one limiting case of a double-base propellant, 
there is not necessarily reason to suppose the double-base pyrolysis law to be applicable to NG 
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Table 3. Values of Condensed-Phase Parameters Used in NG Calculations 


Value 

Description 

As —1.8 x 10 3 g/cm 2 -s 

Exponential prefactor in pyrolysis law 

Es - 9935 cal/mole 

Activation energy in pyrolysis law 

ps = 1.59 g/cm 3 

Mass density of liquid NG 

W NG = 227.087 g/mole 

NG molecular weight 

AHf K =-88,600 cal/mole 

NG heat of formation at 298 K 

cf “ = 0.2975 cal/g-K 

Specific heat of liquid NG 


solely on this account. However, both NG and nitrocellulose are nitrate esters whose likely first step 
in decomposition is an NO 2 scission; if this were a rate limiting step, then the decomposition rates 
for NG and nitrocellulose might be expected to be similar. This argument may, in fact, be the 
underlying reason for the apparent universality of the double-base pyrolysis law. In any case, we 
will have to let the results be the ultimate justification for this assumption. 

A major source of ambiguity in applying the new model to NG is the identification and 
quantification of the condensed-phase decomposition products, which are the nascent gas-phase 
species. As suggested previously, this is a situation in which a number of possibilities (consistent 
with a balanced reaction) must be examined and, ideally, arbitrated by broad theoretical reasoning. 
Again, the results must be the ultimate justification for this choice. We have examined the use of 
three sets of condensed-phase decomposition products for the NG case. The first is motivated by the 
only other elementary reaction treatment of NG of which we are aware. Hatch [11] assumed, 
without attempting justification, that NG decomposed by the path, as indicated in Table 4. 


Table 4. NG Condensed-Phase Decomposition Paths Considered 


Overall NG Decomposition Path 

Reaction Path Label 

NG (C 3 H 5 N 3 O 9 ) □ 3 NO 2 + 2 CH 2 0 + HCO 

Hatch [11] 

NG □ 2 NO 2 + HONO + 2 CH 2 O + CO 

Levy [391 

NG □ 3 HONO + 2 HCO + CO 

MSM4 
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We have not found a way to rationalize this path theoretically, but deem it worthwhile to 
consider it anyway based on its previous use. A second path, given in Table 4, was hypothesized by 
J. B. Levy [39] based on theoretical considerations. The steps in NG decomposition proposed by 
Levy are illustrated in Figure 9. A third set (MSM4) of decomposition products is also given in 
Table 4 and is used to compute burning rates. No theoretical justification based on a detailed 
sequence of paths is offered for this set, but it does lead to good agreement with experimental rates 
(as will be shown). One of our aims in considering a number of condensed-phase reaction paths is to 
gauge the sensitivity of the computed burning rate to the choice of path. 

Figure 10 shows the results of the burning rate calculations using each of the above sets of 
condensed-phase decomposition products at an initial temperature of298 K. The figure also shows 
experimental burning rates measured by Andreev [40, 41] and Andreev et al. [42]. The MSM4 
decomposition product set leads to very good agreement with the experimental burning rates over a 
wide range of pressure. It is also of interest to note that the burning rate can vary by as much as an 
order of magnitude among the possible decomposition product sets that we have considered. This 
fi ndin g could be considerably important to tailoring the propellant burning rate for specific purposes 
if a way could be found to actively influence the decomposition pathways in the condensed phase. 

In all of the calculations presented here involving NG (including those of the next section), the 
windward differencing option and suppression of the thermal diffusion and multicomponent 
transport options were selected. The numerical grid options were such that the burning rates 
presented are numerically accurate to about 2% or better. 

Since this is the most extensive, detailed chemical modeling of NG combustion ever performed, 
a detailed analysis of the gas phase flame for one condition is presented. The case chosen is for pure 
NG at 10 atm at an initial temperature of 298 K, with the assumed MSM4 product set 
(3HONO+2HCO+CX)). This is the product set which yields the best agreement between predicted 
burning rates and experiment (see Figure 10). Pathway diagrams, shown in Figure 11, have been 
constructed for this case using a postprocessor code, PREAD, written at ARL for the PREMIX code. 
To obtain these diagrams, the rates were first integrated over distance from the propellant surface to 
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Figure 9. Mechanism for Condensed-Phase Decomposition of NG Proposed by Levy [39]. 

0.030 cm. This region, which we term the primary reaction zone, includes all of the near surface 
flame up to the leading edge of the dark zone, as determined by plateaus in temperature and key 
species (e.g. NO and concentration profiles [not shown]). Using the integrated rates gives a concise 
global picture of the chemistry in this region. The flame is quite structured in this region, however, 
having five separate peaks in the heat-release profile. A complete discussion of this structure is 
outside the scope of report, but some remarks about the structure are contained in the following 
discussion. The chemistry of the dark zone of NG, the region of low gradients in species 
concentration and temperature between the primary flame and the secondary flame (visible flame), is 
similar to that for double-base propellants. That case has been discussed elsewhere [37,38]. 
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Linear Burning Rate (cm/s) 



Pressure (MPa) 


Figure 10. Computed NG Burning Rates at 298 K for Different Decomposition Paths Using a 
New Semi-Empirical Model Compared With Experimentally Measured Rates of 
Andreev [40,41] and Andreev et al. [42]. 
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(a) Carbon-Chemistry Dominant Pathways. 



HNO 


(b) Nitrogen-Chemistry Dominant Pathways. 

Figure 11. Identification of Dominant Chemical Pathways for NG at 10 atm and 298 KUsing 
the MSM4 Product Set. Numbers in Parentheses Are the Relative Rates of Each 
Reaction Value of 100 in Relative Units, Which Equals 8.817 x 10‘ 3 mol/cm 2 -s. 
Line Thickness Is Proportional to the Total, Relative, and Integrated Rates for 
Each Path. 
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The conversion of HCO to other species is very rapid, which is not surprising considering its 
radical nature. This conversion takes place very close to the surface, occurring within the first 
0.001 cm. As shown in Figure 11a, most of the HCO is converted to CO by unimolecular 
decomposition or by reaction with NO, H, or NO 2 . However, a portion of the HCO is temporarily 
converted to CH 2 O, primarily by reacting with HNO. The CH 2 O, being much less reactive than 
HCO, survives further into the primary flame. Later in the primary reaction zone, the CH 2 O is 
converted back into HCO by reacting with OH and H; the HCO formed in this way never builds to 
an appreciable concentration, but is instead converted to CO or CO 2 first stage products. 

The nitrogen chemistry is considered in Figure 1 lb. As with HCO, the conversion of HONO to 
other species begins rapidly near the surface because of its direct reaction with a portion of the HCO, 
which forms NO 2 and CH 2 O. However, only a small portion of the available HONO reacts prior to 
the depletion of the HCO. Instead, most of the HONO conversion is delayed, occurring mainly 
between 0.001 and 0.015 cm. A portion of the HONO is converted directly to NO via unimolecular 
decomposition. However, much of it is first converted to NO 2 via reaction with OH, H, or HCO. 
Curiously, the conversion of NO 2 to NO takes place mainly after the burnout of the HONO. 
Ordinarily, in combustion involving NO 2 as a major species, it is found that NO 2 converting to NO 
by reacting with H atoms is so fast that it occurs at a very early point, and it controls the radical pool 
concentrations. But here, the HONO+M = NO+OH+M reaction is the major radical source, the 
HONO concentration in the region is large, and the HONO+OH rate constant is also large. This 
leads to HONO+M=NO+OH+M and HONO+OH = H 2 0+N0 2 as the major nitrogen pathways close 
to the surface, rather than NO 2 +H = NO+OH. The NO 2 formed early in the flame from HONO 
reacts higher in the flame after the HONO is gone. N0 2 is converted to NO primarily by an H atom 
reaction, but there are also important contributions from HNO, HCO and CO reactions. NO, the 
only nitrogenous species surviving into the dark zone, maybe viewed as the nitrogenous “product” 
of the first stage combustion. The concentration of HNO is small, except very close to the surface. 
It appears that although their concentrations are small, HNO and NO play important roles in a 
catalytic cycle for the near-surface HCO conversion. 
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More information concerning the importance of the various reactions to the predicted burning 
rate may be obtained from a sensitivity analysis. The sensitivity coefficient of the temperature at a 
given spatial grid point to a particular reaction rate is defined as 


s,-A 

T 


dT 


(14) 


where A, is the A-factor of the ith reaction, T is the temperature at the point of interest, and T m is the 
maximum temperature observed in the entire solution. The second grid point in the solution is 
chosen as the point of interest because this one has special significance. The burning rate is strongly 
controlled by the amount of heat feedback from the gas-phase reactions to the propellant. The heat 
feedback rate is in turn primarily controlled by the temperature gradient at the propellant surface. 
Thus, the magnitude of rate constants for reactions that have a high temperature sensitivity at the 
second grid point in the solution must have a strong influence on the computed burning rate. In fact, 
to a good approximation, the sensitivity coefficient of the heat feedback is proportional to the 
sensitivity coefficient of the second grid point. Table 5 shows the first-order temperature 
sensitivities for the second gas-phase grid point of the solution (and the heat feedback) for reactions 
with relative values above three. For the coefficient of the second grid point, a relative value of 100 
equals an absolute value of7.769 x 10’ 4 . A positive temperature sensitivity coefficient indicates that 
an increase in the rate coefficient of the reaction will result in an increase in temperature at point 
two; such an increase is expected to increase the burning rate. Note that the largest sensitivity is for 
HCO+M = H+CO+M. This is not surprising since the reaction not only converts HCO to a final 
product, but produces the highly reactive H atom. Sensitivities of most of the other reactions can 
also be rationalized by their effects on the radical pool or formation of final products with release of 
heat. The only one whose high, positive sensitivity has been difficult to rationalize is the second one 
on the list, H 2 +NO 2 = HONO+H. This reaction is reversed throughout the entire calculational 
domain (e.g.. Figure lib). Its reversal is not surprising since HONO is initially in the MSM4 
decomposition product set. Ordinarily, since this reversed reaction converts a highly reactive H 
atom to comparatively unreactive species, one would expect the corresponding sensitivity would be 
negative. It seems likely that the explanation for the positive sensitivity is the formation of NO 2 , 


25 




Table 5. Relative Sensitivity Coefficients of the Temperature Gradient at the Surface to the 
Most Sensitive Reactions (in Rank Order) for NG, With MSM4 Product Set at 
10 atm and 298-K Initial Temperature. Positive Values Indicate Reactions That 
Increase the Temperature Gradient When the Rate Coefficient for the Reaction Is 
Increased. 


Reaction Number (see Appendix) 

Relative Sensitivity Coefficient 

33 HCO + M^H + CO + M 

100 

139 H^NO^HONO + H 

46.2 

4 NO + OH + M^HONO + M 

29.4 

34 H + HCO^H 2 + CO 

-26.4 

178 HCO + NCh^H + CCh + NO 

23.6 

175 HCO + HNO ^ CH 2 0 + NO 

21.3 

177 HCO + NO^ CO + HONO 

-21.3 

176 CH 2 0 + N0 2 ^HC0 + H0N0 

19.6 

128 H + HNO^H 2 + NO 

-15.9 

147 HCO + NO ^HNO + CO 

-4.97 

157 H + CH 2 O^HCO + H 2 

-3.94 

64 N02 + H ^ NO + OH 

3.17 


which in turn leads to NO and HNO formation close to the surface. The NO and HNO take part in a 
catalytic cycle for rapid HCO conversion to CO and CH 2 0 in this region (see Figure 11, both parts). 
This conversion occurs via reactions between HCO+NO = HNO+CO and HCO+HNO = CH 2 O+NO, 
which release considerable heat near the surface. 


Rate constants of most of the highly sensitive reactions shown in Table 5 are well established. 
Three exceptions are R175 (HCO+HNO=CH 2 0+NO), R178 (HC0+N0 2 =H+C0 2 +N0), and R177 
(HCO+N02 = CO+HONO). The rate constant expression for R175 has especially large error limits. 
The estimate is attributed to Tsang and Herron [43], who obtained it based on comparison to similar 
reactions. Rate constant expressions for R177 and R178 were obtained from Lin et al. [44]. The 
expressions were computed based upon ab initio and TST-RRKM theory. This approach should 
yield reasonable estimates, but experimentation is desirable. 
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7. Effect of Chemical Additives on the Burning Rate 


It has always been hoped that theoretical modeling might some day contribute to a solution for 
the problem of the effects of chemical additives on the burning rate of propellants. However, only 
with the relatively recent advent of chemically specific modeling with elementary reactions was 
there any real prospect for realizing these hopes. In this section we demonstrate that theoretical 
guidance useful to the propellant formulator, namely the effect of additives on the burning rate, is 
becoming feasible. 

As a first demonstration of the effects of additives on propellant flames, we took the converged 
values of burning rate (0.394 cm/s) and surface temperature (628 K), obtained in a calculation of NG 
at 10 atm using the Levy decomposition product set and the DB11 reaction mechanism, as our 
starting conditions for a series of steady premixed flames with several different additives. A 10% 
mole fraction of NH 3 was then added to the Levy product set (reducing the original product mole 
fractions proportionately). The flame structure was then computed and compared to the pure NG 
case—the result is shown in Figure 12. The two-stage flame structure is known to occur because of 
the reduction of HONO and NO 2 to NO in the primary reaction zone (next to the surface) and the 
subsequent, slower reduction of NO to N 2 in the secondary reaction zone (at the end of the dark 
zone). The secondary reaction zone is also known as the secondary flame or the visible flame. In 
that figure, the secondary gas flame for the pure NG case stands off from the surface by about a 
centimeter. (Notice that the distance scale is logarithmic.) This is typical of the dark zone length of 
double-base propellants. When the NH 3 is added, the dark zone collapses by about a factor of five. 
Also, the heat feedback increases by 31%. Normally, when the heat feedback increases, one 
can expect the burning rate to increase as well. The effect NH 3 has on the dark zone length may 
explain why M30 (29% nitrocellulose, 22% nitroglycerine, 47% nitroguanidine, 2% stabilizer) bums 
with no apparent dark zone, unlike any other gun propellant. A major ingredient in M30 is 
nitroguanidine, which may be expected to supply NH 2 upon decomposition. The chemical rationale 
for this effect will be discussed later in this report. A case with 10% N 2 added is shown in Figure 12 
for comparison. The N 2 has its expected diluent effect, lowering the heat feedback by 28%, but has 
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Distance from Surface (cm) 

Figure 12. The Effect of Two Additives, NH 3 and N 2 , on the Flame Structure of NG at 10 atm 
and 298 K Using Levy’s [39] Decomposition Product Set Surface Temperature 
and Mass Flux Are Fixed at 628 K and 0.626 g/cm 2 -s, Respectively, for All Cases. 

relatively little effect on the dark zone length. Liau [15] performed a similar study of the effect of 
various additives on the dark zone of RDX. 

A more systematic approach is needed to compare the effects of one additive to another in a 
quantitative sense. Adopting the propellant formulation strategy of using enough additive to bring 
the mixture to a zero oxygen balance affords a rational method of determining the amount of each 
additive appropriate for comparison purposes. The oxygen balance is defined as that amount of 
oxygen one must add or subtract to have all oxygen appear in either H 2 0 or C0 2 . NG has a positive 
balance of 3.5% (i.e., it has an oxygen surplus). Thus, we compute that one third of a mole of NH 3 
must be added to effect a neutral oxygen balance for each mole of NG. For comparison, we examine 
two other potential additives, CH 2 0 and H 2 . We are, of course, limited here to those fuel molecules 
that are already in our reaction mechanism. 
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Additional assumptions are required in order to model the effects of condensed-phase mixtures 
of ingredients. The mass density of the mixture is computed by the method of additive partial molar 
volumes, such as 


Pm = W -»/X X M • 


(15) 


where W avg is the average molecular weight, X, is the mole fraction of ingredient i, and V,- is the 
molar volume of ingredient i. This quantity is important in calculating the linear burning rate from 
the mass burning rate. Secondly, we assume that the starting mixture enthalpy is given by the 
weighted sum of the ingredient enthalpies. This is a reasonable approximation to make, but it 
ignores any enthalpy of mixing or solution contributions; there is also some ambiguity as to what 
state is best to use for the additive enthalpy. Finally, we must, of course, assume that the pyrolysis 
law is unchanged as a result of these additives. This assumption is reasonable in view of the small 
amount of each additive. 

The results of the additive computations are given in Table 6 for the Levy decomposition product 
set and Table 7 for the MSM4 decomposition product set. The enthalpies of all of the additives in 
Table 6 and Table 7 at 298 K are taken as that for the gas phase. It should be noted that this 
assumption could have important consequences for the computed burning rate. For example, if the 
enthalpy for NH 3 liquid at 298 K is used, the linear burning rate increases by only 12%. This smaller 
increase reflects the additional energy requirement of the heat of vaporization. Clearly, the 
calculation of the unreacted mixture enthalpy may require a more sophisticated theoretical treatment. 
In interpreting the small percentage differences shown in these tables, bear in mind that the 
numerical accuracy of these calculations is no worse than 2 %; thus, some of the smaller effects 
found may not be significant. 

Through sensitivity and rate analysis, we were able to identify the chemical mechanism for the 
burning-rate enhancement of NG by NH 3 . Essentially, NH 2 reduces NO to final product N 2 by the 
following two most important paths, 
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Table 6. Effect of Different Chemical Additives on the Burning Rate of NG at 10 atm and 
298 K Assuming the Levy [39] Decomposition Product Set 


Additive 

Weight % 

Linear Burning Rate 
(% increase) 

Mass Burning Rate 
(% increase) 

Heat Feedback 
(% increase) 

nh 3 

2.4 

19 

15 

19 

ch 2 o 

3.2 

1 

-1 

-5 

h 2 

0.44 

1 

-8 

-6 


Table 7. Effect of Different Chemical Additives on the Burning Rate of NG at 10 atm and 
298 K Assuming the MSM4 Decomposition Product Set 


Additive 

Weight % 

Linear Burning Rate 
(% increase) 

Mass Burning Rate 
(% increase) 

Heat Feedback 
(% increase) 

nh 3 

2.4 

16 

12 

14 

ch 2 o 

3.2 

7 

5 

2 

h 2 

0.44 

2 

-7 

-6 


NH 2 + NO -» NNH + OH 

n 2 + h 2 o. 

The NNH thus formed rapidly decomposes to N 2 + H. The reaction therefore leads to chain 
branching, which maybe important to the burning rate enhancement, in addition to the increased NO 
conversion rate to N 2 . The NH 2 is formed from NH 3 , predominately by the following reaction: 


NH 3 + OH -* NH 2 + H 2 0. 


This mechanism reinforces our speculation that M30 forms no dark zone due to the presence of an 
amino group on the nitroguanidine. 
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8. Speculations on Practical Burning-Rate Modifiers 


The mechanism of NH 3 action invites a further discussion of the M30 case. With 
NH 3 -modified NG, we found both a collapse in the dark-zone length and a significant boost in the 
burning rate. M30 has no dark zone, consistent with our NH 3 -modified NG case, but it bums no 
faster than a single-base propellant and slower than a double-base propellant, a fact which appears to 
be inconsistent with our model calculation. This apparent inconsistency might be explained as 
follows. It is likely that the large percentage of nitroguanidine (NQ) present in M30 is responsible 
for the lack of burning-rate enhancement since it acts to lower the flame temperature of the 
propellant. Hence, the cooling effect might outweigh the rate acceleration effect for such major 
proportions of NQ. We would expect that if a small amount of NQ, say 2-5%, was added to either a 
single or double-base propellant, then the burning rate would be increased and the dark-zone length 
diminished. Of course, some of the benefit of the NH 2 from NQ could be hindered by the energy 
required to either break down its parent molecule or to get it into the gas phase, as could be seen in 
our model calculation’s sensitivity to the NH 3 starting enthalpy state. Nonetheless, here is a 
concrete, theoretically inspired idea that could easily be tested. Thus, though the model needs to be 
further refined and expanded, it can already provide insights of potential worth to the propellant 
formulator. 

The propellant formulator may find it expedient to tailor decreases in the burning rate as well as 
increases. For example, it has been found [45] that for layered propellants, at least a three-to-one 
ratio in the rates of the two layer materials must be attained in order to achieve desired performance 
enhancements. We have conducted preliminary calculations [46] showing HNCO as an additive that 
both decreases the heat feedback (and therefore presumably the burning rate) and decreases the 
visible flame standoff. As NH 2 proved to be the active molecule in increasing the heat feedback 
resulting from the addition of NH3 to NG, NCO appears to be the active molecule for the HNCO 
additive. One factor leading to the reduced heat feedback may be that the radical NCO combines 
with NO and NO2 near the surface to form relatively stable molecules such as N2, N2O, CO, and 
CO2, thereby diminishing the rate of radical growth. Reduction in the visible flame standoff (dark 
zone length) may be due to the unimolecular decomposition in the dark zone of the N 2 O which is 
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formed in the near surface region via the NCO reactions. This decomposition, N 2 O -» N 2 + O, 
increases the radical growth rate in the dark zone, causing acceleration of the reactions leading to the 
visible flame. In fact, a decrease in the visible flame standoff may be a desirable effect in that it 
might reduce gun ignition delays, delays attributable to the slowness of the reactions leading to the 
substantial heat release of the visible flame [47]. Of course, to realize the benefits of these 
provocative ideas, much work is yet to be done both in searching for suitable active molecules and in 
finding practical substances that can release these active molecules. Now, the way is clear for a new 
and productive partnership between the theoreticians and the formulation chemists. 

9. Conclusions 


Models of energetic material combustion have increased enormously in their sophistication and 
complexity in recent years. Treatments of the gas phase involving dozens of species and hundreds of 
elementary reactions are now common. On the other hand, chemically specific descriptions of the 
condensed phase have not kept pace, largely because of the lack of definitive experiments and theory 
relating to the reaction paths there. All existing treatments of the condensed phase consider several 
overall reactions at most, and even these suffer from uncertainties in the identities and concentrations 
of the reactants and products. In the face of these intractabilities, we have developed a new semi- 
empirical model, one which requires as input a single overall reaction representing the conversion of 
unreacted energetic material to the nascent gas-phase species and a pyrolysis law expressing the 
relationship of surface temperature to burning rate. This approach, where the pyrolysis law is 
universal over some class of materials, makes the most sense; the existence of such a relation has 
been previously verified for double-base propellants. We have demonstrated that the simplest form 
of the pyrolysis law also well describes surface regression due to evaporative and reactive 
mechanisms. The capability of this model in describing the burning rate as a function of pressure 
and even species profiles (where available) was illustrated for frozen ozone, RDX, and 
nitroglycerine. 

Key uncertainties in applying this model to energetic materials include the identity and 
concentrations of the products of the condensed-phase decomposition. By examining a number of 
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possible product sets, we found that the computed burning rate could vary by as much as an order of 
magnitude among these sets. This finding suggests that it might be possible to influence the burning 
rate of a material by stimulating decomposition along normally less favored paths. 

One of the long-held hopes for chemically specific combustion modeling is the capability to 
compute the effects of chemical additives on the burning rate, in order to rationalize the performance 
tailoring of propellants. We believe that progress is rapidly being made towards fulfilling that 
promise. As a demonstration of this budding capability, we determined the relative effect of three 
additives on the burning rate of nitroglycerine where the amounts of each additive were computed to 
achieve a neutral oxygen balance for each mixture. 
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N02+0=NO+02 


3.900E+12 

0.00 

-238.0 

66. 

NCO+H=NH+CO 


5.4 00E+13 

0.00 

0.0 

67. 

NCO+0==NO+CO 


4.52 OE+13 

0.00 

0.0 

68. 

NCO+N=N2 +CO 


2.000E+13 

0.00 

0.0 

69. 

NCO+OH==NO+CO+H 


2.000E+13 

0.00 

7500.0 

70 . 

NCO+M=N+CO+M 

N20 

H20 

N2 

C02 

Enhanced by 
Enhanced by 
Enhanced by 
Enhanced by 

1.140E+23 

5.000E+00 

5.000E+00 

1. OOOE+OO 

1.500E+00 

-1.95 

59930.0 

71. 

NC0+N0=N20+C0 


8.80 0E+17 

-1.78 

790.0 

72 . 

NC0+N0=C02+N2 


1.130E+18 

-1.78 

790.0 

73 . 

NC0+N02=C02+N20 


1.950E+13 

-0.26 

-620.0 

74 . 

NCO+N02=CO+NO+NO 


1.770E+12 

-0.26 

-620.0 

75. 

NH+02=HN0+0 


4.610E+05 

2.00 

6500.0 

76. 

NH+02 =NO+OH 


1.2 80E+06 

1.50 

100.0 

77. 

NH+N0=N20+H 


3.500E+14 

-0.46 

16.1 

78. 

NH+NO=N2 +OH 


2.160E+13 

-0.23 

0.0 

79. 

N20+H=N2 +OH 
Declared duplicate 

reaction... 

2.53 0E+10 

0.00 

4550.0 

80. 

N20+H=N2 +OH 
Declared duplicate 

reaction... 

2.23 OE+14 

0.00 

16750.0 

81. 

NNH+0=N20+H 


1.4 00E+14 

-0.40 

477.0 

82 . 

NNH+0=NO+NH 


3.3 00E+14 

-0.23 

-1013.0 

83 . 

N20+0=N2 +02 


3.654E+12 

0.00 

15900.0 

84 . 

N20+0=N0+N0 


9.985E+13 

0.00 

28040.0 

85. 

H+HNO=NH+OH 


3.OOOE+14 

0.00 

18000.0 

86. 

NH+0H=N+H20 


5.000E+11 

0.50 

2000.0 

87.. 

NH+N=N2+H 


3.000E+13 

0.00 

0.0 

88. 

N+H2-NH+H 


1.600E+14 

0.00 

25140.0 

89. 

HN0+H=NH2+0 


3.500E+15 

-0.30 

28200.0 

90. 

NH2 +0=NH+0H 


6.750E+12 

0.00 

0.0 

91. 

NH2+0H=NH+H20 


4.000E+06 

2.00 

1000.0 

92 . 

NH2 +H=NH+H2 


4.000E+13 

0.00 

3650.0 

93 . 

NH2 +NH=N2H2 +H - 


1.500E+15 

-0.50 

0.0 

94 . 

NH2+N=N2 +H+H 


7.20 0E+13 

0.00 

0.0 

95. 

NH2 +02 =HNO+OH 


4.50 OE+12 

0.00 

25000.0 

96. 

NH2+NH2=N2H2+H2 


5.000E+11 

0.00 

0.0 

97. 

NH2 +NH2 =NH+NHS 


5.000E+13 

0.00 

10000.0 

98. 

NH2+NH2=N2H3+H 


1.790E+13 

-0.35 

11320.0 

99. 

NH2+NH2 +M=N2H4+M 


2.980E+47 

-9.44 

9680.0 

100. 

NH2+N02=N20+H20 


2.840E+18 

-2.20 

0.0 

101. 

NH+N02=N20+OH 


1.000E+13 

0.00 

0.0 

102 . 

N2H4 +H=N2H3 +H2 


1.OOOE+12 

0.50 

2000.0 

103 . 

N2H4 +OH=N2H3 +H20 


3.000E+10 

0.68 

1290.0 

104. 

N2H4 +0=N2H3 +OH 


2.000E+13 

0.00 

1000.0 

105. 

N2H3 =N2H2+H 


1.200E+13 

0.00 

58000.0 

106. 

N2H3 +H=N2H2 +H2 


1.000E+12 

0.50 

2000.0 

107. 

N2H3 +OH=N2H2 +H20 


3.000E+10 

0.68 

1290.0 

108. 

N2H3 +0=N2H2+0H 


2.000E+13 

0.00 

1000.0 

109. 

N2H2+M=NNH+H+M 

H20 

02 

N2 

H2 

Enhanced by 
Enhanced by 
Enhanced by 
Enhanced by 

5.000E+16 

1.500E+01 

2.000E+00 

2.000E+00 

2.000E+00 

0.00 

50000.0 

110. 

N2H2 +H=NNH+H2 


5.000E+13 

0.00 

1000.0 

111. 

N2H2 +0=NH2 +NO 


1.000E+13 

0.00 

0.0 
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112. N2H2 +0=NNH+0H 

113. N2H2 +0H=NNH+H20 

114. N2H2 +NH=NNH+NH2 

115. N2H2 +NH2 =NH3 +NNH 

116. NH2 +NO=NNH+OH 

117. NH2+N0=N2+H20 

118. NH3+0H=NH2+H20 

119. NH3 +H=NH2 +H2 

120. NH3 +0=NH2 +OH 

121. NH3 +M=NH2 +H+M 

122. NNH+NO=N2+HNO 

123. NNH+H=N2+H2 

124. NNH+OH=N2 +H2 O 

125. NNH+NH2 =N2 +NH3 

126. NNH+NH=N2+NH2 

127. HN0+0H=N0+H20 

128. H+HNO=H2 +NO 

129. HNO+NH2-NH3+NO 

130. N+NO=N2+0 

131. 0+NO=N+02 

132. NO+H=N+OH 

133. HN0+HN0==N20+H20 

134. N20+N0=N2+N02 

135. NO+NO+NO=N2O+N02 

136. HOCO+M=OH+CO+M 

137. C0+N02 =N0+C02 
13 8. CH+N02 =HCO+NO 

139. H2 +N02 =HONO+H 

140. NNH=N2+H 

Declared duplicate reaction... 

141. NNH+M=N2+H+M 

Declared duplicate reaction... 

142. HN0+N0+N0=HNN0+N02 

143. HNNO+NO=NNH+N02 

144. HNNO+NO==N2 +H0N0 

145. HNNO+M=H+N20+M- 

146. HNNO+M=N2+OH+M 

147. HCO+NO=HNO+CO 

148. 0+CH20<=>0H+HC0 

149. 0+CH20H<=>0H+CH20 

150. 0+CH30< = >OH+CH2 O 

151. 0+CH30H<=>0H+CH20H 

152. 0+CH30H<=>0H+CH30 

153. 02+CH2 0< = >HO2 +HCO 

154. H+HCO (+M)< = >CH20(+M) 


Low pressure limit: 

0.13500E+25 

TROE centering: 

0.78240E+00 

H2 

Enhanced by 

H20 

Enhanced by 

CO 

Enhanced by 

C02 

Enhanced by 


2.000E+13 

0.00 

1000.0 

1.000E+13 

0.00 

1000.0 

1.000E+13 

0.00 

1000.0 

1.000E+13 

0.00 

1000.0 

9.3 OOE+11 

0.00 

0.0 

2.OOOE+2 0 

-2.60 

924.0 

2.040E+06 

2.04 

566.0 

5.420E+05 

2.40 

9917.0 

9.400E+06 

1.94 

6460.0 

2.200E+16 

0.00 

93470.0 

2.000E+13 

0.00 

0.0 

1.000E+14 

0.00 

0.0 

5.000E+13 

0.00 

0.0 

5.000E+13 

0.00 

0.0 

5.000E+13 

0.00 

0.0 

1.295E+07 

1.88 

-958.0 

4.460E+11 

0.72 

655.0 

2.OOOE+13 

0.00 

1000.0 

3.270E+12 

0.30 

0.0 

3.800E+09 

1.00 

41375.0 

1.700E+14 

0.00 

48800.0 

3.63 0E-03 

3.98 

1190.0 

4.290E+13 

0.00 

47130.0 

1.070E+10 

0.00 

26800.0 

2.190E+23 

-1.89 

35270.0 

9.040E+13 

0.00 

33780.0 

1.010E+14 

0.00 

0.0 

3.210E+12 

0.00 

28810.0 

3.000E+08 

0.00 

0.0 

1.OOOE+13 

0.50 

3060.0 

1.7 00E+11 

0.00 

2100.0 

3.200E+12 

0.00 

270.0 

2.600E+11 

0.00 

810.0 

2.200E+15 

0.00 

21600.0 

1.000E+15 

0.00 

25600.0 

7.230E+12 

0.00 

0.0 

3.900E+13 

0.00 

3540.0 

1.OOOE+13 

0.00 

0.0 

1.OOOE+13 

0.00 

0.0 

3.880E+05 

2.50 

3100.0 

1.3 00E+05 

2.50 

5000.0 

1.000E+14 

0.00 

40000.0 

1.090E+12 

0.48 

-260.0 


-0.25700E+01 0.14250E+04 

0.27100E+O3 0.2755 0E+04 0.65700E+04 

2 . OOOE+OO 
6. OOOE+OO 
1.500E+00 
2 .OOOE+OO 


155. H+CH20 (+M) < = >CH20H (+M) 


5.400E+11 0.45 3600.0 
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Low pressure limit: 

0.12 700E+33 

-0.48200E+01 0.653 0 0E + 04 



TROE centering: 

0.71870E+00 

0.103 0 0E + 03 0.12 910E + 04 0 

. 41600E+04 


H2 

Enhanced by 

2.000E+00 




H20 

Enhanced by 

6.OOOE+OO 




CO 

Enhanced by 

1.500E+00 




CO 2 

Enhanced by 

2.000E+00 



156 . 

H+CH20(+M)<=>CH30(+M) 

5.400E+11 

0.45 

2600.0 


Low pressure limit: 

0.22000E+31 

-0.48000E+01 0.55600E+04 



TROE centering: 

0.75800E+00 

0.94000E+02 0.1555 0E+04 0 

. 42000E+04 


H2 

Enhanced by 

2.OOOE+OO 




H20 

Enhanced by 

6.000E+00 




CO 

Enhanced by 

1.500E+0 0 




C02 

Enhanced by 

2.000E+00 



157. 

H+CH20<=>HC0+H2 


2.300E+10 

1.05 

3275.0 

158. 

H+CH20H(+M)<=>CH30H(+M) 

1.800E+13 

0.00 

0.0 


Low pressure limit: 

0.3 0000E+32 

-0.48000E+01 0.33000E+04 



TROE centering: 

0.76790E+00 

0.33800E+03 0.18120E+04 0 

. 50810E+04 


H2 

Enhanced by 

2.000E+00 




H20 

Enhanced by 

6.000E+00 




CO 

Enhanced by 

1.5 00E+00 




C02 

Enhanced by 

2.000E+00 



159. 

H+CH20H<=>H2+CH20 


2.000E+13 

0.00 

0.0 

160. 

H+CH30(+M)<=>CH30H(+M) 

5.000E+13 

0.00 

0.0 


Low pressure limit: 

0.86000E+29 

-0.40000E+01 0.30250E+04 



TROE centering: 

0.8902 0E+0 0 

0.14400E+03 0.28380E+04 0 

. 45569E+05 


H2 

Enhanced by 

2.000E+00 




H20 

Enhanced by 

6.000E+00 




CO 

Enhanced by 

1.500E+00 




C02 

Enhanced by 

2.OOOE+OO 



161. 

H+CH30<=>H+CH20H 


3.400E+06 

1.60 

0.0 

162 . 

H+CH30<=>H2+CH20 


2.000E+13 

0.00 

0.0 

163 . 

H+CH3 OH< = > CH2 OH+H2 


1.70 OE+07 

2.10 

4870.0 

164 . 

H+CH30H<=>CH30+H2 


4.200E+06 

2.10 

4870.0 

165. 

H2+CO(+M)<=>CH20(+M) 


4.30OE+07 

1.50 

79600.0 


Low pressure limit: 

0.50700E+28 

”0.34200E+01 0.84350E+05 



TROE centerings 

0.93200E+00 

0.1970 0E+03 0.1540' 

0E+04 0 

. 103 00E + 05 


H2 

Enhanced by 

2.OOOE+OO 




H20 

Enhanced by 

6.000E+00 




CO 

Enhanced by 

1.5 00E+00 




C02 

Enhanced by 

2.OOOE+OO 



166. 

0H+CH20<=>HC0+H20 


3.430E+09 

1.18 

-447.0 

167 . 

0H+CH20H<=>H20+CH20 


5.000E+12 

0.00 

0.0 

168. 

0H+CH30<=>H20+CH20 


5.000E+12 

0.00 

0.0 

169. 

0H+CH30H<=>CH20H+H20 

1.440E+06 

2.00 

-840.0 

170 . 

OH+CH30H< = >CH3 0+H2 0 


6.300E+06 

2.00 

1500.0 

171. 

H02+CH20<=>HC0+H202 


1.000E+12 

0.00 

8000.0 

172 . 

CH+H20<=>H+CH20 


1.713E+13 

0.00 

-755.0 

173 . 

CH20H+02<=>H02+CH2O 


1.800E+13 

0.00 

900.0 

174 . 

CH30+02 < = >H02+CH20 


4.28 0E-13 

7.60 

-3530.0 

175. 

HC0+HN0==CH20+N0 


6.000E+11 

0.00 

2000.0 

176. 

CH2 0+N02 =HCO+HONO 


8.020E+02 

2.77 

13730.0 

177. 

HC0+N02 =CO+HONO 


1.240E+23 

-3.29 

2355.0 

178 . 

HC0+N02=H+C02+NO 


8.390E+15 

-0.75 

1930.0 
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Note: Units for the rate parameters are centimeter, second, and mole, and for E, cal/mole. For 
reactions followed by three numerical parameters, the rate-coefficient expression is 

L 

k = AT exp(-E/RT). For reactions which appear twice with the phrase “declared duplicate 
reaction...,” the rate coefficient is computed as the sum of the two three-parameter expressions. For 
reactions involving a generalized collider species, M, collider efficiencies other than 1.0 are 
specified. For reactions involving pressure-dependent rate expressions (those with a collider species 
specified as [+M]), three types of expression are used. If “T&H VALUE” occurs in the output, the 
Tsang and Herron form was used (as described in Tsang and Herron* 1 ), with constants ao and ai (if 
the latter is used) appearing, respectively, on the same line. A version of CHEMKIN, modified at 
ARL, was used for this computation. If “TROE centering:” occurs, the TROE form was used with 
the appropriate parameters specified on that line. If neither of these is mentioned, the Lindemann 
form was assumed. Descriptions of the TROE and Lindemann expressions may be found in the 
CHEMKIN manual. 2 


* Note that the log expressions used in this source are for base 10 (W. Tsang, private communication). 

1 Tsang, W., and J. T. Herron. “Chemical Kinetic Data Base for Propellant Combustion I Reactions Involving NO, N0 2 , 
HNO, HNO 2 , HCN, and N 2 0.” Journal of Physical and Chemical Reference Data , vol. 20, p. 609-663,1991. 

2 Kee, R. J., F. M. Rupley, and J. A. Miller. “Chemkin-II: A Fortran Chemical Kinetics Package for the Analysis of 
Gas-Phase Chemical Kinetics.” Sandia National Laboratories Report SAND89-8009, September 1989. 
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